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Abstract. Two key reactions of hydrostatic nuclear burning in stars have recently been revised by new 
experimental data - the '"'Nfp, and 3 a reaction rates. We investigate how much the new rates in¬ 
fluence the evolution of low-mass, metal-poor and metal-free, stars and of an intermediate-mass star of 
solar-type composition. We concentrate on phases of helium ignition or thermally unstable helium burn¬ 
ing. Our global result is that the new 3 a rate has no significant influence on such stars, but that there is a 
noticeable, though small effect of the new y)*^0 rate, in particular on the core helium flash and the 

blue loop during core helium burning in the intermediate-mass star. 

Key words. Stars: interiors - Stars: evolution - Nuclear reactions 


Send offprint requests to: A. Weiss; 
(e-mail: weiss@mpa-garching.mpg.de) 
Correspondence to: A. Weiss 



2 


Weiss et al.: Stellar models with revised reaction rates 


1. Introduction 


The triple-alpha reaction is one of the key nuclear reactions for the synthesis of the elements in stars and is 
also the main energy source during helium burning. The rea ction rate is d ominated by resonances, the best 
known being the one at 7.65 MeV, theoretically predicted b}? Hovlellll9.54ll. but there is considerable interest 

bo et al.l t 


in determining all resonances with high precision. Recentlv. lFvnbo etaf l20n.5ll reported new measurements 
concerning resonances of with 3 a-particles obtained from ^^C-decay experiments. In particular, they 
found a dominant re sonance at ^ 11 MeV, while they did not confirm another at 9.1 MeV reported previously 
( lAngulo et alJll999ll. The new re action rate (called the “ISOL” rate in the following) deviates from the one 
published by Angulo et alJ (Il999l NACRE): at temperatures between 2.5 ■ 10^ and 10* K it is between 7 and 
20% lower, but below 2.5 ■ 10^ K up to 50% higher. The largest change occurs for T > 3 ■ 10® K, where the 
new rate is increasingly lower by up to one order of magnitude (Fig.^. 

In the present paper we consider three cases of stellar evolution where a modified 3 a rate might in¬ 
fluence the models: the core helium flash in low-mass metal-poor stars, the core and shell helium burning 
in intermediate-mass stars of solar-like composition, and helium burning in metal-free low-mass stars. This 
mass range is interesting, because of the high temperature sensitivity during core and shell flashes, where 
instabilities amplify even small temperature variations. 

Coincidentally, another key reaction has been re-determined recently by the LUNA collaboration, the 
7)^^0 bottleneck reaction of the CNO-cyclej|Formicolajetaj200^. The influence on globular cluster 
age determinations has already been discussed bv llmbrianreraiT ( )2004l) . and we will briefly comment on this. 
Our main interest, however, in this bottleneck reaction is again how stars of low and intermediate mass might 
be influenced in their evolution. 

The outline of this paper is as follows: After describing briefly our stellar evolution code and the specific 
reaction rates for the 3 a process and the ^'^N reaction in Sect. 2, we present the results for the three stellar 
evolution cases described above in Sect. 3, followed by the conclusion in the final section. 
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Fig. 1. 3 a reaction rates relative to the N ACRE analytic one dAng iilo et al. |l999l) . Shown are the tabulated 
NACRE, CFHZ85 llCaughlan et al.l 198.51) and ISOL rates (iFvnbo et alJl2005l) . The dotted line is the absolute 
value of the reference NACRE rate (scale on right axis). 
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2. Stellar evolution program 

2.1. Basic properties 


All calculations were done with the Garching stellar evolution code as described bv lWeiss & Schlattll ll200r)ll . 
The program was modified in only minor aspects since th en. To summarize b riefly the main physical ingre- 
dients, the code incorpora tes the OPAL equation of state ilRogers et a QIIM and the OPAL opacit y tables 
dlglesias & Rogerslll996l) supplemented by the molecular opacities of Afexander^F^gusonl lll994l) . In the 
calculations presented here, the equation of state (EOS) of Irwin (see ICass'isretan]2003ll has been used, 
which is based on the OPAL EOS. Microscopic diffusion of all elements is implemented and has been in¬ 
cluded in the calculations of Sect. l3.2l Convection is treated according to standard mixing le ngth theory wit h 
the Schwarzschild criteri on for stability. Mass l oss is included according to Reimers’ formula ilReimersll97-5h . 
in the generalized form dlben & Renzinil[r983ll containing a free scaling parameter rj, which we will specify 
for each case. 

Nuclear burning is taken into account by a network which treats hydrogen and helium burning separately, 
unless both protons and a-particles are present and temperatures are high enough for helium processing. In 
such cases (see Sect. the w hole network and any mixing process are treated simultaneously. Eor details 
on this see lScTdattl et al1ll200ll) . 


2.2. and 3 a rate 

The standard implementation of the first rate' is according to the recommendation in lAdelberger et al ] (ll998i 
Table VI; hereafter referred to as Adel98), with the astrophysical 5(0)-factor being 3.5 keV b. NAC RE 
dAngulo et alJll99^ gives S(0) - 3.2 + 0.8 keV b. The new LUNA value dEormicola et al.l 120041) is 
1.7 + 0.1 (stat) + 0.2 (sys) keV b; S'(0) and S”{0) have been left unchanged at the Adelberger values and 
th e new rate therefore d iffers from the old one by a constant factor. This is in agreement with the procedure 
of llmbriani et al.ld2004l) . 

The standard implementation of the 3 a rate in our program is the analytic form of ICaughlan et alJ d 19851 
CEHZ85). Since lEvnbo et alJ d2005l) used the analytic fit to the NACRE rate for reference, we implemented 
the same rate for the comparisons. The new ISOL rate was available to us in tabulated form (Chr. Diget, 
private communication). Eig.^shows all 3 a rates compared with the NACRE analytic rate. N ote that in the 
tempe rature range of interest to us the new rate lies within the uncertainty of the NACRE rate dAngulo et al.l 
I 1999 I) and that the analytic fit of the NACRE rate devia tes from the tabulated values by up to 5%. The ISOL 
rate agrees quite well with that of lCaughlan et al.ldl98-5ll over an extended range below temperatures of several 
billion degree. The critical temperature above which helium burning is taken into account was set to 5 • 10'' K; 
using still lower values has no influence on the models. 

3. Sample calculations 

3.1. Core heiium fiash in a Pop. II star 

In this section we present the results of our sample calculations and compare them with reference calcu¬ 
lations that use the NACRE and Adel98 rates for the 3 a and *''N(/7,y)*^0 reaction rates, respectively. We 
begin with the onset of helium burning in stars of low mass and degenerate helium cores, the so-called 
core helium flash. We performed calculations for two masses, 0.8 Mq and LOMq, and three metallicities, 
Z = 0.001, 0.0001, 0.00001 and initial helium content Y - 0.25. No mass loss was assumed. The calculations 
were started from the zero-age main sequence (ZAMS) and followed through the flash until the star settles on 
the horizontal branch (HB) and then up to the early AGB. 

Updating the 3 a rate, the changes are very small but as would be expected from a slightly lower rate. 
The tip of the RGB is almost unchanged, the luminosity being 2% higher for the new rate. The core mass 
increases very little, between 0.002 and 0.003 Mq depending on mass and metallicity. The age of the models 
also increases by a mere 0.2 Myr at most. We show the resulting tracks in the Hertzsprung-Russell Diagram 
(HRD) in Eig.|2]for the lowest metallicity and mass, and list a few quantities in Tabled Note that even the 
secondary flashes during the settling to the HB are almost indistinguishable. The agreement between the two 
cases is as close for all other cases. 

We then varied the *"'N(p, 7)'^0 rate {3 a - NACRE). The effect turned out to be much larger, and can 
be seen in Eig. |2l (lower panel). We notice differences already around the turn-off, and in particular during 

‘ T he term “rate” is to be understood as **{crv), that is as the Maxwellian-averaged reaction rate (see lAngiilo et al.l 
I 1999 I) in units of cm^ mol^' s“' for n reaction partners. 
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Fig. 2. Evolution of a stellar model of 0.8 Mq and Z = 10 from ZAMS to HB. In the upper panel we varied 
the 3 a- and in the lower one the 7)'^0-rate. The insets show the approach to the HB in greater detail. 


the approach to the HB; the effective temperature of the zero-age HB model increases by up to 20% for the 

LUNA rate. _ _ 

At this point it is worthwhile to clarify the result found bv llmbriani et alJ(l2004ll . They showed that due to 
the improved '^N(p, y)*^0 rate the ages of globular clusters increase by up to 1 Gyr. This is not simply due to 
the slower bottleneck reaction leading to longer timescale for hydrogen fusion. In Fig.|3we display for one 
example (1 Mq, Z = 0.001) the evolution of luminosity and effective temperature during main-sequence and 
red-giant-branch evolution. Evidently, the luminosity evolution on the main sequence is hardly affected by the 
nuclear rate. Energy generation takes place within a core that has a temperature only marginally higher but is 
more extended as the result of a shallower temperature profile. So, while at the center hydrogen is burning at 
a slower pace (same T, lower reaction rate), more mass is involved in the nuclear fusion. The time it takes to 
deplete hydrogen in the core is therefore somewhat higher (in our case h- 190 Myr), but the exhausted core is 
larger (-1-0.013 Mg). Due to that larger core mass, towards the end of the main sequence the central temperature 
is higher by approximately 5% for the LUNA rate. Another consequence of the different core evolution is a 
modified temperature gradient, and thus at the end of the main sequence evolves differently (Fig. |3, 
influencing the turn-off luminosity which is defined as the value of L at the bluest point along the MS track. 
To be more quantitative, the turn-off in the case of the Adelberger rate is at f = 4.38 Gyr, log Tefr = 3.852, 
and logL/Lo = 0.554. For the revised LUNA rate the hottest point is reached later, at log Tefr = 3.854, when 
logL/Lg = 0.593 and t — 4.62 Gyr, an increase of 240 Myr. Note, however, that at the turn-off age of the 
Adelberger case, — 3.854 and logL/Lg = 0.558 for LUNA, i.e., almost identical. Thus, the change in 
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Table 1. Selected properties of one of our Pop. II models (M = 0.8 Mq; Z = 10 at critical stages of its 
evolution from ZAMS to ZAHB evolved with various combinations of reaction rates as shown in Fig. ElMne 
is the helium core mass in solar unit at the helium flash. 




Adel98 

LUNA 

3a: 


NACRE 

NACRE 

ISOL 


logL/Lo 

0.437 

0.467 

0.467 

Turn off 

log Teff 

3.8360 

3.8368 

3.8368 


Age [Myr] 

11130 

11273 

11273 


log L/Lq 

3.285 

3.237 

3.246 

RGB-tip 

log Teff 
Afue 

3.6449 

0.4990 

3.6479 

0.5052 

3.6474 

0.5074 


Age [Myr] 

12638 

12699 

12699 

ZAHB 

log LI Lq 

1.744 

1.715 

1.722 

log Teff 

3.8614 

3.9324 

3.9344 


the turn-off is due to the modihed interior structure, which influences the morphology of the evolutionary 
track. We have also investigated a few other cases and find that at even lower metallicity (Z = 0.0001) the age 
differences decrease to 1-2%. 
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Fig. 3. Comparison of luminosity and effective temperature as functions of time (in Myr) for a star with 1 Mq 
and Z = 0.001, calculated with both the Adelberger (dotted line) and LUNA (dashed line) '"^N(p, 7)^^0 rate. 
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3.2. Core helium flash in a Pop. Ill star 

A particularly interesting variant of the core helium flash is that found i n met al-free (Pop. Ill) stars. Details 
on this subject were provided bv lWeiss et all hoool) and lSchlattl et al.l 1I2OOII1 . We have calculated the evo- 
luti on of a star of M = 1 M n. F, = 0.23 (Z = 0). We varied both rates as before, the reference case using 
the lAdelberger et ^ ( Il998l) and NACRE 3 a rates. Some properties of the models are listed 

in Table |2] We also investigated the case “LUNA + CFHZ85”, but do not discuss it further, because it is 
indistinguishable from the “LUNA + ISOL” case. 

The first particularity of Pop. Ill evolution is a loop in the HRD after the main sequence, which is due 
to spu rious carbon production in the hot core and the ignition of the CNO-cycle (for details see lWeiss et alJ 
200nl) . The presence of this loop depends on the stellar mass; it is present in our reference case, but vanishes 
in all cases with the new (lower) LUNA rate for ^'*N(/7,y)'^0. As the reason for the disappearance of the 
loop we identify the fact that the lower CNO-luminosity prevents the creation of a convective core, which 
otherwise leads to an increase of the hydrogen abundance in the core and thus to a luminosity enhancement. 
We note that the creation of the first in-situ carbon nuclei is not influenced by the choice of the 3 a rate. 


Table 2. Selected properties of Pop. Ill models evolved with various combinations of reaction rates. We list 
luminosity and helium core mass at the onset of the first and second core helium flash, as well as the maximum 
helium luminosity (logL|^“/Lo) and the mass at which the maximum temperature is reached (Mj^^IMq) at 
the time of logL^^/Lo. In addition, the surface abundances (mass fractions) after the flash-induced mixing 
event is given for key elements. 


'^N-t-p 


Adel98 

LUNA 

30- 


NACRE 

NACRE 

ISOL 

MS 

CN-loop 

yes 

no 

no 


log(L/Lo) 

2.324 

2.302 

2.311 

P‘ flash 

log(LS“/Lo) 

10.08 

10.01 

10.03 

Mne/Mo 

0.4784 

0.4755 

0.4780 



0.1849 

0.1818 

0.1863 


H 

0.5288 

0.5369 

0.5366 


He 

0.4582 

0.4510 

0.4510 

dredge up 


0.00417 

0.00292 

0.00302 

surface 

'^C 

0.00111 

0.00080 

0.00083 

abdcs. 

14IVJ 

0.00772 

0.00837 

0.00861 


isjsi 

2.6 ■ 10'^ 

1.1 • 10-'' 

1.1 • 10'^ 


160 

3.9 • 10'® 

3.4 • 10'® 

4.0 • 10'® 


log(L/Lo) 

3.532 

3.487 

3.498 

2 "'* flash 

log(LSr/Lo) 

MhcIMq 

9.592 

0.4585 

9.764 

0.4615 

9.770 

0.4638 



0.0789 

0.0665 

0.0729 


The second particular event is the (first) core helium flash that happens at a much lower luminosity in 
Pop. Ill stars than in those with Z > 10'®, due to the higher core temperature. As the H-shell temperatures are 
higher as well, the so-called entropy barrier is lower and a mixi ng event between the hydrogen-rich envelope 
and the hot carbon-enriched helium layers is possible (see also lFuiimoto et a] hood . As a consequence, the 
envelope will be enriched drastically in CNO-products and helium (see Table|2l. As in the case of ordinary 
Pop. II stars, the new 3 a rate has only a very moderate influence on the flash: it starts at slightly higher 
luminosity and core mass, but the increase is not significant. The peak helium luminosity is unchanged, 
within the computational uncertainties. 

In contrast, the new ^'*N(/7,y)'®0 rate has a much stronger influence. The flash luminosity is lower due 
to the higher shell temperatures and the core mass is higher. This is probably a consequence of the fact that 
the core mass-luminosity relation depends on the temperature exponent of the shell hydrogen burning. For 
CNO-burning, it is L cx Mj, and for pp-burning oc M^. Very metal-poor stars have a very small contribution 
from the CNO-cycle only and therefore follow closely the latter relation. A lower *"^N(p, 'y)*®0 rate leads to a 
smaller temperature exponent, and thus the core mass is slightly higher at given luminosity. 
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Fig. 4. Evolution in the HRD for the Pop. I model with M - 5 Mq, showing the influence of changing the 
*'*N(p,y)'^0 reaction rate. The two cases are indicated in the figure. 


The result of the flash-induced dredge-up shows larger (relative) variations (Table|2ji, but with no apparent 
systematics, except that in the table the dredge-up appears generally to be stronger for the older rates, while 
for the latest rate more CN-burning could take place. 

After the star has settled again on the giant branch, it resumes its evolution since core helium burning 
has been extinguished as a result of the first flash and mixing. A second, more moderate flash appears at 
a standard RGB tip luminosity (logL/L q 3.5). Table |2 demonstrates again that the largest changes are 
due to the new LUNA *^N(/7,7)'^0 rate, leading to lower ignition luminosity but higher helium core mass. 
The largest variation found in Table |2is therefore between the reference case (col. 3) and that updating the 
'"^N(p, y)'^© rate only (col. 4). 

3.3. Helium burning and thermal pulses in a Pop. I star 

As a third case we have investigated that of a typical solar-metallicity intermediate-mass star evolving into and 
through the Asymptotic Giant Branch (AGB) phase. The composition of our model was X = 0.695, Y - 0.285 
and Z = 0.02. Mass loss according to Reimers was taken into account with the scaling parameter j] being 0.4 
on the RGB and 0.5 on the AGB. 

In Fig. 0]we show the evolution of this model in the HRD for the two *'*N(p, y)*^0 reaction rates. Using 
the updated LUNA rate the main-sequence evolution takes place at higher temperatures, implying a more 
compact structure of the model. This is similar to the Pop. II model (Sect. 3.1) but more pronounced. As the 
evolution of intermediate mass stars is very sensitive to the internal composition profile, the consequences for 
the core helium burning phase are quite large; the blue loop gets significantly shorter. 

The luminosity variations during the thermal pulses (TP) along the AGB are shown in Fig. |3 The TPs 
start earlier by about 1 million years for the LUNA rate, which is due to a shorter main-sequence lifetime. We 
have therefore shifted the time axis for the Adel98 case to coincide at the pre-flash luminosity maximum. The 
pulse behaviour is similar, although the peak luminosity and peak helium luminosity (not shown) are higher, 
and the interpulse duration longer for the new LUN A rate. This is consisten t with earlier findings concerning 
the influence of H-burning rates on pulse behaviour jPespain & Scaloll976l) . and in particular with the recent 
investigation bv iHerwig & Austir] ll2004h . who found stronger flashes for lower *^N(p, 'y)*^0 rates in a 2 Mq 
model. These differences become smaller with increasing pulse number. We followed 15 (LUNA) respectively 
8 (Adel98) pulses and then stopped the calculations. 

Using the new LUNA *"^N(/7,7)^^0 reaction rate, we then compared models calculated with the three 
variants of the 3 a rate. Again, the blue loops, being most sensitive to structure variations, show the clearest 
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Fig. 5. Thermal pulses on the AGB for two different reaction rates. Shown is the total lumi¬ 

nosity as function of time (in million years). To display better the differences the time axis for the Adel98 
case has been shifted by -997 kyr. 


reaction, and get even shorter for the new ISOL rate. Concerning the thermal pulses, they are extremely similar 
to each other, and if shifted by -80 kyr for the CFHZ85 and -67 kyr for the new ISOL 3 a rate relative to the 
NACRE case they completely agree. 

Finally, we mention the chemical composition of all cases investigated in terms of central carbon and 
oxygen abundance after core helium burning. The carbon abundance ranges from 0.200 to 0.223 with the 
highest value reached for the older reaction rates. Consequently the core most rich in oxygen is obtained 
for the LUNA *'*N(p, y)*^0 and new 3 a rates combined. However, in spite of noticeable variations, they are 
minor compared with the uncertain ties still present due to the '^C(a, 7)'®0 rate. As an example, we quote 
here the result of iBono et alJ ll2000l their Tab. 6) that the time spent in the blue loop by a 5 M© star changes 
by up to 2.12 Myr (33%), the mass of the CO-core by 0.26 Mq (5%), and the C/O ratio by a factor o f 3 when 
varying the rate within a factor of 2.3, which is comparable to the rate uncertainty dKiinz et a itooA 


4. Conclusions 

We have investigated the influence of new rates for the *'^N(p, y)*^© and 3 a reactions on the evolution of low- 
and intermediate-mass stars, considering cases in which helium burning proceeds under thermally unstable 
conditions. In all cases we find negligible changes in the evolution and also in the interior evolution due to 
the new 3 a rate. The RGB tip brightness is slightly increased due to the lower rate at very low temperatures 
corresponding to the earliest phases of helium ignition. The largest effect shows up in the blue loops during 
core helium burning of the 5 Mq star, emphasizing the sensitivity of these loops to details of the interior 
structure. 

The '^N(p, y)^^0 rate has a definitely stronger influence. It prolongs the duration of central hydrogen 
burning, increases the turn-off temperature and thus indirectly the turn-off location, it leads to a disappearance 
of the CN flash in the Pop. Ill post-main-sequence star, strongly reduces the blue loop during core helium 
burning of a 5 Mq star, and also influences the thermal pulses. Interestingly, these are, apart from a very minor 
shift in time, almost identical in the case of varying the 3 a rate. 

We conclude that the new 3 a rate has no influence on the evolution of low and intermediate-mass stars, 
and that the effect of the LUNA y)'^0 is tiny, but noticeable. Further investigations into its effect on 

more massive stars as well as of the effect of the 3 a rate on massive stars are indicated. If the differences in 
the various rate determinations we have used are representative of the experimental errors, then both rates are 
no longer a source of uncertainty for stellar modeling. 
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